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TOWARDS THE PROOF OF YOSHIDA’S CONJECTURE 


Jim JAHN AND JIRINA JAHNOVA 

Abstract. Yoshida’s Conjecture formulated by H. Yoshida in 1989 states that in equipped with 
the canonical symplectic form dp A dq, the Hamiltonian flow corresponding to the Hamiltonian 

^ N N 

H with qo = qn+i = 0, 

i=l i=0 

where iV > 3 is odd and fc > 4 is even, has no global complex meromorphic first integral functionally 
independent of H. For TV = 3 and TV = 5 with fc > 4 arbitrary even number, the result was proved true 
by Maciejewski, Przybylska and Yoshida in 2012 by means of differential Galois theory. However, the 
question whether Yoshida’s conjecture is true in general, remained open. In this paper we give a proof 
that this conjecture is in fact true for infinitely many values of TV using the results of R. D. Costin which 
are based on the so-called poly-Painleve method devised by M. Kruskal. 


1. Introduction 

General Hamiltonian dynamical systems are traditionally of great interest for both mathematicians 
and physicists simply because they describe the evolution of many physical systems like a planetary 
system, electron in an electromagnetic held etc. Since the number of independent hrst integrals de¬ 
termines the dimension of the space that the trajectories hll and therefore provides some information 
about the possible onset of chaos in the phase space, it is no surprise that one of the basic problems 
concerning these systems is to decide whether a given system is integrable or not in the sense of Liouville 
theorem (cf. [1]), which means to decide whether there are sufficiently many functionally independent 
hrst integrals in involution that are of a given class or not. In the realm of integrability theory, it is 
a standard trick of the trade already since the days of S. Kowalevski to consider complexihed systems 
instead of the original ones (see e.g. UDIEIDSIEIEIIIIIEIIEIEIES]) and virtually the same route 
is taken also in this paper - here, we consider Hamiltonian systems in the complex domain (the phase 
space is a subspace of and the time variable is complex). In the 80’s, Ziglin’s theory of complex 
nonintegrability of these complexihed systems, in which the necessary conditions of integrability are 
expressed in terms of certain monodromy matrices of the solutions of the corresponding variational 
equations along a suitably chosen particular solution of the original Hamiltonian system, started to de¬ 
velop [23], [25]. This theory was further extended by H. Yoshida [2l], [22]. Building on Ziglin’s theory, 
Morales-Ruiz and Ramis were able to obtain further results using the diherential Galois theory that 
relate the integrability of Hamiltonian systems to the solvability of their variational equations around 
a particular solution (see [13], [H], [15] and HZ])- The necessary conditions of partial integrability of 
Hamiltonian systems by means of diherential Galois theory have been recently studied by Maciejewski, 
Przybylska and Yoshida HD], HB. 
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Although the theory of (non)integrabilhy is developing, there are still systems that seem not to be 
amenable to analysis using this theory. One of these systems is the N-degrees-of-freedom Hamiltonian 
system 


dq 

dt 


dH 


dt 


dH 


q= (gi,...,g7v) e C , p = (pi, ... ,pjv) e C 


N 


whose dynamics being given by the Hamiltonian 


1 


N 


H = + vw,fc(q), 


( 1 ) 


( 2 ) 


2=1 

where the potential is the following homogeneous function 

N 

^iv,fc(q) = - Qi+if, with go = Qn+i = 0, (3) 

i=0 

and the time t is considered to be a complex variable. 

To the best of our knowledge, this system was first considered by H. Yoshida in [22] where it was 
shown it does not admit a single additional first integral for k = 4 and N = 3,5 except the Hamiltonian 
itself. In this paper Yoshida formulated the following conjecture 


Yoshida’s Conjecture. For an arbitrary odd N > 3, and for arbitrary even k > 4, the Hamiltonian 
system with N degrees of freedom, given by the Hamiltonian ^ with the potential ^ does not admit 
an additional global complex meromorphic first integral. 


Yoshida himself actually proved that the conjecture is in fact a corollary to another conjecture of him 
which he also stated in [22] in the following form: 


For an arbitrary odd N >3 and even k > 4, the numbers Ai,..., Aat given by 

= + 1<]<N. (4) 

where Xj := 2{k — 1) sin^((7rj)/(2(Y + 1))), are independent over Q. 


However, in 2007 K. Yoshimura proved that this latter conjecture is not true when he found a 
counterexample to the case N = 5, k = 16 (for details see [23l |TT]), so that the question whether 
Yoshida’s Conjecture is true or not remained open. In 2012, Maciejewski, Przybylska and Yoshida m 
proved that Yoshida’s Conjecture is true for N = 3,5 and arbitrary even k > 4 using the results obtained 
in n] that are based on an application of differential Galois theory to variational equations along a 
particular solution, and they strongly conjectured that Yoshida’s Conjecture is true for an arbitrary 
odd N > 3 and even k > 4. Unfortunately, the approach through differential Galois theory seems to be 
virtually impossible to be applied to systems that contain parameters, so that the problem of finding 
the general proof of Yoshida’s Conjecture is still open. 

The principal aim of this paper is to give a proof of Yoshida’s Conjecture for infinitely many values 
of N, more precisely for those N > 3 such that N = 1 mod 6 or A = 3 mod 6 with an arbitrary even 
k > 4. To this end we employ results due to R. D. Costin which are of slightly different nature than the 
results discussed above. They are based on the ideas of Martin D. Kruskal’s poly-Painleve method [6], 
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0. i. 0 which, although at times somewhat heuristic in spirit, hnds its completely satisfactory and 
rigorous counterpart in the work of Costin. These results first appeared in |2] and their generalization 
was later given in the Ph.D. thesis of Costin (see HI)- 

The present paper is organized in the following way: the main ideas and results pertaining to the 
proof are described in Section [2l To streamline the proof, we briefly recall several basic facts about 
Chebyshev polynomials in Section [3l these will be later used in the proof. The very proof of Yoshida’s 
Conjecture for the special values of N cited above occupies the rest of the paper and is to be found in 
Section HI 

From now on, the term “(non)integrability” stands for “complex (non)integrability”, and “meromor- 
phic function” means “complex meromorphic function”. 

2. The poly-Painleve test: a review of results due to R. Costin 


In this section, we give a short account of the results related to polynomial homogeneous systems of 
second-order ordinary differential equations introduced in [ 3 ] that are a generalization of similar results 
obtained in p] in the case of the homogeneous Henon-Heiles system and at the same time they are 
elaboration of the poly-Painleve method devised by M. D. Kruskal. 

Consider the system of second-order differential equations 


gm = Rn(q), m = l,...,N, q = (gi,. .. ,gAr), f e C (5) 

with Pm homogeneous polynomials of degree k — 1. Particular solutions qm{t) of ([5]) that are of the form 
Qmit) = am4>it), m = 1,..., N, am G C, satisfy the conditions 

(f) (f) , Otm • • • 1 •Tv) (6) 

for every m = 1,..., Y. Let a = (ai ,... ^a^) be a nonzero solution of ([ 6 ]) such that ciAr 7 ^ 0. If we 
introduce new variables Mi, ..., mw-i and Q, a small parameter e by the formulas 

(y. 

qm{t) = —Q{t) + eUm{t), m = l,...,Y-l, 
an 

qN(,t) = Q{t) 


and if we Taylor expand the polynomials Pm, then using homogeneity, the system ([5]) takes the form 

ii = Q’^~'^a]^’^~'~'^Mpu +eR{u,Q,e), (7) 

Q = Q’^-^a]^’^^^PNia)+eRN{u,Q,e) ( 8 ) 


where u = (ui,..., mat-i), i? is a vector of polynomials, Rn is a polynomial and Mp is a matrix, given 
by the formula 


Mp = 



Ti dPjsf , A 

dqj J.. 




(9) 


If the matrix Mp has N — 1 distinct eigenvalues Ai,..., Atv-i, then the reduced system (i.e. the system 
corresponding to equations ([7]), ([ 8 ]) with e = 0 ) can be converted, after integrating once the reduced 
form of equation (IH]) and eliminating time by treating Q as an independent variable and rescaling Q, 
into the decoupled system of Y — 1 generalized Lame equations of the form 


- 1 ) 




m , rv dUf] 




k 

+ -X- 


dx 


-1^ 


x^ ‘^v 

m^’j <-'n 


= 0 , 


( 10 ) 
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for m = 1,..., — 1. 

In |1], it is proved that r functionally independent first integrals Fi,..., of the system ([S]) dehned 
on a domain D = F x C such that 1) the projection Fg^^ of T) on the ^Tv-coordinate contains 

closed paths around the roots of the polynomial — iCk) (as a polynomial in q) for some 
constant F G C and 2) all F^’s are meromorphic along the linear manifold qm = ^Qn, Qm = ^Qn, 
m = — 1, give rise to r independent hrst integrals of the system ([TU]) that are moreover 

holomorphic on a domain F whose projection on the a:-coordinate contains closed paths around all the 
Fth roots of unity. 

At this point, employing the results concerning such hrst integrals of the system flTUD . Costin was 
able to prove the following theorem: 

Theorem 1 (R. Costin, |1]). Suppose the matrix Mp has N — 1 distinct eigenvalues Ai, ..., Aat-i. Let 
n numbers among the numbers 

{k - 2)2 k m = l,...,N -1 (11) 

16 2 ^ ^ 

be irrational. Then the system ^ has at most 2N — 1 — 2n independent first integrals which are 
meromorphic near the linear manifold q^ = Qm = ^Qn, rn = 1,..., N — 1. 

3. ChEBYSHEV polynomials of first and second KIND: A BRIEF REVIEW 

In this section, we briehy recall several rudimentary facts concerning the Chebyshev polynomials of 
the hrst and second kind. These technical preliminaries will be used in Section 0] where the certain cases 
of Yoshida’s Conjecture are proved. A general reference for this section is na or [19]. If not specihed 
otherwise, every occurrence of the symbols n, m in this section refers to arbitrary nonnegative integers 
n and m. 

One of the possibilities to dehne the Chebyshev polynomial of the first kind Tn{x) is to require that 
it is the unique solution of the recurrence relation 

T„+i(a;) = 2xTn{x) - Tn-i{x) (n > 1), (12) 

with initial conditions To{x) = 1, Ti{x) = x. The Chebyshev polynomials of the hrst kind satisfy the 
following functional equation, called the nesting property. 

Tn {Tm{x)) = Tnm{x). (13) 

The Chebyshev polynomial of the second kind Un{x) may be uniquely dehned as solutions of the recur¬ 
rence relation 

Un+l{x) = 2xUn{x) - Un-l{x) (u > 1) (14) 

with initial conditions Uo{x) = 1, Ui{x) = 2x. The following representation for the Chebyshev polyno¬ 
mial of the second kind is standard: 




, n > 1. 


(15) 
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Next we list a handful of properties that show the mutual connection between the Chebyshev polyno¬ 
mials of the first and second kind and that will turn out to be useful in the sequel: 

Tn{x) = ]^{Un{x)-Un- 2 {x)) for every u > 2, (16) 

n 

Un{x) = 2 Tj{x), if n is odd, (17) 

j odd 
n 

Un{x) = 2 Tj{x) — 1, if n is even. (18) 

j=0 
j even 

Finally, we shall later need the following transformation property between the Chebyshev polynomials 
of the hrst and second kind: 


fomn— i(t) Ura—l iTn{x^^ Un—\{x^ ^ (?7l, 77. ^ 1). 


(19) 


Another important sort of problems, especially when dealing with polynomials, is the question of the 
location and nature of their roots. The Chebyshev polynomials of either kind are a well-explored area 
in this respect, since both Tn{x) and Un{x) have exactly n distinct roots lying in the interval (—1,1) 
that are given by the formulas 


x1 = cos 


Ti{2k — 1) 
2n 


Xu = cos 


kiT 

n -|- 1 


( 20 ) 


for fc = 1 ,..., 77. 


4. Main results 


Consider the Hamiltonian system given by ([2]). It can be shown that the system ([2]) can be equivalently 
written as a second-order system of the following form: 

^ = -VViv,»(q), (21) 

where V stands for the usual gradient operator in cartesian coordinates. We note that — VVAf,A:(q) 
is a homogeneous polynomial of degree k — 1, since is a homogeneous polynomial of degree k in 
q = (gi,..., gTv)- Hence the results from Section [2] are applicable in this case. 

The corresponding point a = (cii,..., ajq) that describes the particular solutions of the form qm{t) = 
for the system fl2ll) is any solution of the algebraic system 


^ —ka\ ^ — k{ai — 02)^ ^ ^ 


( a, \ 

k{ai — 0:2)^”^ — k{a2 — 0:3)^“^ 

= 

0(2 

^ /c(ajv_i - y 


Utv y 


( 22 ) 


Using the ansatz a = (ui,..., uat) = (a, 0, —a, 0,..., (—1)^“a), we obtain one particular solution a 
of the system fl22ll with a = {—l)^^{2k)^^, where (—1)*^ is in general complex-valued, taken with 
principal values of the argument. 
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In this particular case, the {N — 1) x [N — 1) matrix Mp from Q corresponding to P = —Wjsf^k 
takes the following form: 


Mp = -(k - 1) 


/ 

2 

-1 

0 

0 

... 0 



( ° 

0 

0 . 

. 0 

ai 

\ 


-1 

2 

-1 

0 

... 0 



0 

0 

0 . 

. 0 

02 



0 

-1 

2 

-1 

... 0 


1 

0 

0 

0 . 

. 0 

O3 

















0 

0 

-1 

2 



2 ^ ’ 

0 

0 

0 ■ 









-1 






0 

O-N-2 


V 

0 

0 

0 


-1 2 




0 

0 . 

. 0 

ON-1 

/ 


where the numbers a,- read like this: 


(23) 


CLi — 


(—1)*"^^ ^ if f is odd, 


0 


if i is even. 


We have the following 

Lemma 1. The characteristic polynomial XpW of the matrix Mp from 1^) is equal to 

1 


2N- 




(24) 


where x = 2 — 


Proof. Clearly, the characteristic polynomial xp of fho matrix Mp in fl23|) is given by the determinant 
of the matrix 


Mp-XE 


:(*:-l) 


/ 

X 

-1 

0 .. 

. 0 

bi \ 


-1 

X 

-1 .. 

. 0 

b2 


0 

-1 

X 




0 

0 

-1 

. -1 

1 

CO 





X 

bN-2 

V 

0 

0 

0 . 

. -1 

X / 


(25) 


with X = 2 — bi = Oi ioT i = 1,2,N — 3, bM -2 = —2. The polynomial Xp(A) = det {Mp — XE) 
can be computed by the Laplace expansion along the {N — l)st column: 


Xp(A) 


^(fc -1)~-‘ x;(-i)"-'+‘5»(-i)"-'-‘A-. 

k=l 

N-l 


( 26 ) 
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where /,• denotes the determinant of the j x j sub-matrix 


/ 

X 

-1 

0 

0 

... 0 

\ 


-1 

X 

-1 

0 

... 0 



0 

-1 

X 

-1 

... 0 



0 

0 

-1 

X 








-1 



0 

0 

0 


— 1 X 



(27) 


for j = 0,1,..., N — 2 with /q = 1. This is a Toeplitz tridiagonal matrix whose determinant is the 
so called continuant with both the minor diagonals consisting solely of (—l)’s and with the diagonal 
consisting of x’s. It is a well known fact (cf. [H]) that the continuant satisfies the following three term 
recurrence relation 

/j+i = xfj - (28) 

which in our case moreover satisfies the initial conditions /o = 1 and /i = x. This is exactly the recur¬ 
rence relation satished by [/j(|)’s from flTT)) . Thus Xp(A) is equal to the following linear combination 
of Chebyshev polynomials of the second kind: 

N-l 


Xp(A) — 



(29) 


2fc+iV-l .. 

2 U‘ 


k=0 


2k 


2k+N-l (X 


X 


77 ) ~ ‘^ Un -3 77 ) + xUn -2 77 


(30) 


the second one results from plugging 

N-5 


the hrst equality follows upon substituting Uk (|) for fk into 
the corresponding bkS into fl29l) (the resulting sum with the upper index is to be interpreted as 
void for N = 3) and the last equality is just a consequence of the definition flTTl) . At this point we split 
the proof in two parts depending on whether the number of summands of the sum from 
(the first case) or odd (the second case). 

First case. This case is equivalent to the fact that N = 3 mod 4. Then the sum from 
rewritten in the form 


IS even 


can be 


iV-1 

2 




2fc+iV-l 

-^U2k 


k=0 


= Un-1 


— Un-3 


+ fAjV-5 


- . . . + f/2 


-t/n 


= 2 


--If 


X 


+ Tjv-5 t: + • • ■ + ^2 t: 


X 


(31) 


where we used flTbl) in the second equality. Since N is odd, thus — 1 is even, we have N — 1 = 2r for 
certain r G N, r odd, and we can use flT^ to recast the expression fl^ in the following way: 


X 


= Ur 


2 

T2 


X 


X 


2 Tjv-i - +T,v-5 77 +---+^2 77 =2 T, T2 - +T,_2 T2 - +... + T 1 T 2 - 


X 


X 


X 


= Un^ 
2 


T 2 
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Second case. This case is in turn equivalent to the fact that = 1 mod 4. The sum from fl5UD now 
takes the form 


N-1 

2 




k=0 


= 2 T, 


'N-l 


+ ^A-5 


+ . . . + T4 


+ f/o 


(32) 


where again the property ffTBD was employed. Since N is odd, thus — 1 is even, so that N — 1 = 2r 
for certain r G N, r > 2, r even, and using the property ffT3l) we obtain 


2 T; 


'-N-l 


+ 'Tn-s 

X 


+... + T 4 


X 


+ Uo 


= 2 Tr T 2 - +Tr.2 T 2 - +...+T 2 T 2 - 


X 


X 


+ Uo - 


X 


= Ur Tr 


= Ur T, 


2To ( T, 


X 


+ I + U 0 


X 


2 + 1 + 1 


= Un-1 To hr 


X 


where in the second and third equality we respectively used the property 
and Uq{x) = 1 for all x. 


and the fact that To(a;) = 1 

□ 


Remark 1. Note that the transformation property flT^ implies that, upon writing 


X 


X 


Un{- =U,^_,{- =U^{T2{- U, I- =U^{T2{- 


X 


X 


X 


X, 


the characteristic polynomial Xp(A) can be written as 


XpW = 




-U. 


N 


1 - 


A 


k-1 


2 ^- 1 ' ' X 2^{k-l-X) 

To be able to apply Theorem[T], we have to show that the matrix Mp corresponding to P = — VVAr,A:(q) 
has N — 1 distinct eigenvalues. This task is achieved in the following result. 


Proposition 1. The matrix Mp from 
and every even k > 4. 


has precisely N — 1 distinct eigenvalues for every odd N > 3 


Proof. The claim follows from Lemma [U and the last remark since Un{x/2) has precisely N distinct 
roots, whence we can infer that the polynomial 

Mil 

X 

has precisely N — 1 distinct roots. □ 

The next matter of interest that we heavily rely on later is the number of (ir)rational eigenvalues. 
This is treated in 


Proposition 2. The only possible rational eigenvalues of the matrix Mp from are equal to —1) 
or \{k — 1). Moreover, all the eigenvalues belong to the interval {0,2k — 2) for every odd N > 3 and 
every even k > 4. 
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Proof. First we note that the characteristic polynomial Xp(A) of Mp can be written in the form 

N-1 

fc =0 ^ ^ 

In particular, if we normalize the polynomial xp a polynomial in x (which does not affect the roots 
anyway), the leading coefficient of xp as a polynomial of the variable x is equal to 1 and the constant 
term of the same polynomial (of the variable x) is (—1)"^which is always nonzero. Therefore, 
according to the standard result on the rational roots of a polynomial with integer coefficients, we 
conclude that the only rational roots of the polynomial xp ^ are the divisors of (which is always 
integer, since N is odd). In fact, from Remark [1] we can see that all the roots of xp ia x are also 
the roots of the polynomial Un{x/2) and therefore they must belong to the interval (—2, 2). Hence the 
only rational roots of xp as a polynomial of the variable x are the numbers —1,0,1. Using Remark 
[U again, we can easily see that 0 is in fact not a root of xp as a polynomial in x, since x = 0 is the 
only root of f/i(x/ 2 ) = x and therefore, all roots of Un{x/2) being simple, it cannot be the root of 
U n -1 (T-?(x/2)). Since the number x = 2 — is rational if and only if the number A is rational, we can 
infer that the only rational roots of the polynomial Xp(A) are those numbers A satisfying the relations 
±1 = 2 — (2A)/(/c — 1), which means that A = (fc — l)/2 or A = 3{k — l )/2 as claimed. 

The fact that all eigenvalues A belong to the interval (0, 2k — 2 ) clearly follows as x = 2 — (2A )/{k — 1) 
and X G (—2, 2) in view of the considerations made above. □ 


Proposition 3. Suppose N > 3 is odd. Then 

(1) For an arbitrary N such that N = 1 or N = 3 mod 6 and for every even k > A, all the numbers 
Pm in m, m = 1,..., A^ —1, corresponding to the Hamiltonian system given by the Hamiltonian 

are irrational numbers. 

(2) For an arbitrary N such that N = 5 mod 6 and for every even k > A, at least N — 3 numbers 
among the numbers Vm in m, m = l,...,iV — 1, corresponding to the Hamiltonian system 
given by the Hamiltonian are irrational numbers. 


Proof. First note that all the numbers Vm are in fact real (this follows from the previous proposition). 
Then, clearly, if some of the eigenvalues Am of the matrix Mp is irrational for certain m = 1,..., — 1, 

then the corresponding number Vm is also irrational. This means that the only possibility for Vm to be 
rational is when the corresponding number Am is rational. This situation occurs, due to Proposition [21 
only if A = (fc —1)/2 or A = 3{k — \)/2. The fact that these two values of A are roots of the characteristic 
polynomial xp is equivalent to the fact that 


Un^ 

2 




= 0 . 


This last condition is in turn satished only if there is an integer j 


1 

2 


cos 


2j7r 
N +l 


1 ,..., such that 


according to formula ([I]). This happens precisely when N is such that A^ + 1 is divisible by 3 which 
means, since N is odd, that N = 5 mod 6. 








10 


Jim JAHN AND JIRINA JAHNOVA 


□ 

In view of Proposition [3l we have thus proved the following theorem: 

Theorem 2. For an arbitrary N > 3 such that N = 1 or N = 3 mod 6 and for arbitrary even /c > 4, the 
Hamiltonian system with N degrees of freedom, given by the Hamiltonian ^ with the potential ^ does 
not admit an additional first integral meromorphic in a complex neighbourhood of the linear manifold 
Qm = CLmgN, Pm = CLmPN, m = 1,..., N — 1, where Om = (—1) 2 for m odd, and Om = 0 for m even. 
For an arbitrary N such that N = 5 mod 6, the Hamiltonian system with N degrees of freedom, given 
by the Hamiltonian ^ with the potential has at most four additional first integrals meromorphic 
in a complex neighbourhood of the linear manifold = amQN, Pm = (^mPN,fn = 1,... — 1, where 

o-m = for m odd, and = 0 for m even. 

Proof. The claim follows from Proposition [3] and Theorem [H □ 

This last theorem shows that Yoshida’s Conjecture is true for every N such that iV = 1 mod 6 
or N = 3 mod 6 with /c > 4 an arbitrary even number. The case N = 5 mod 6 still leaves certain 
possibility for four additional global meromorphic hrst integrals to exist. However, it is clear that the 
only z/^’s that are possibly rational (and for which the criterion of Costin’s is not decisive) are the 
numbers 

\/5/c^ — 8k + 4 \/13fc^ — 16k + 4 

z/ = -, or u = -. 

4 4 

But these numbers are clearly rational if and only if the corresponding expressions under the square-root 
sign, — 8fc -|- 4 and 13fc^ — 16k + 4, are perfect squares, respectively. Although the task to determine 
precisely those /c’s for which this situation occurs is quite cumbersome (if not impossible), reducing 
our ability to prove Yoshida’s Conjecture in its full generality to a seemingly simple number-theoretic 
problem, the computer-aided numerical experiments that we made suggest that at least among the hrst 
10® fc’s the only even ones, for which the corresponding z/’s are rational, are quite rare (there are in fact 
only 3 of them in each case). Moreover, it turns out that in neither of these cases both z/’s are rational 
simultaneously. This actually seizes the space left for possible additional meromorphic hrst integrals of 
the system in question from 4 to 2 within the range k = 4, ..., 10®. 


k 

V5k'^-8k+4 

^ 4 

Vl3fc^-16Zc+4 

4 

16 

rational 

irrational 

40 

irrational 

rational 

760 

irrational 

rational 

4896 

rational 

irrational 

1576240 

rational 

irrational 

66354520 

irrational 

rational 


Taking into account the results of Maciejewski, Przybylska and Yoshida m , we can state the following 


Corollary 1. Yoshida’s Conjecture holds true: 
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(1) for every N > 3 such that N = 1 mod 6 or N = 3 mod 6 and every even k > 4] 

(2) for N = 5 and for arbitrary even A; > 4; 

(3) for every iV > 11 such that N = 5 mod 6 and every even k = 4,..., 10®, where k ^ 16, 40, 760, 
4896, 1576240, 66354520. 

For every > 11 such that N = 5 mod 6 with k equal to any of the numbers 

16,40, 760,4896,1576240,66354520, 

the Hamiltonian system with N degrees of freedom given by the Hamiltonian ^ with the potential ^ 
admits at most two additional first integrals that are meromorphic in a complex neighbourhood of the 
linear manifold Qm = amqN, Pm = o-mPN, m = 1,..., iV — 1, where Om = (— for m odd, and 
Am = 0 for m even. 

Remark 2. Note that in the case of the system ([2]) there is the following relation between the Yoshida’s 
numbers Aj from (jll) and the Costin’s numbers Vj\ Let Ai,..., \n-i denote the eigenvalues of the matrix 
(l23l) . Then they are equal (according to Remark [T]) to those N — 1 numbers Xj from (j4]) that are different 
from k — 1. Obviously, 



for every j = l,...,iV — 1. Therefore, we can see that instead of independence of Aj’s over Q we can 
require their irrationality as the sufficient condition for non-integrability of the system (|2]), which can 
be sometimes (especially if some parameters are present in the system) easier to decide. 

Remark 3. In general, the fact that a complex Hamiltonian system does not admit an additional 
meromorphic hrst integral (locally or globally) does not imply that the system in question exhibits chaos 
in the real phase space. Therefore, it is interesting to consider a family (indexed by the parameters N 
and k) of real versions of Hamiltonian systems ([2]) and to perform further investigation in this direction. 
One of the tools commonly used in this area is the so-called Poincare section method which rests on 
the study of the discrete dynamical system given by the Poincare map corresponding to a periodic 
trajectory of the system and to a chosen hypersurface in the phase space (the Poincare hypersurface) 
which is transverse to the trajectory. Since it is very difficult to hnd the Poincare map in general, 
some numerical methods are often used in order to hnd the orbit of a given point under this map. 
The orbit can be easily visualized in case of two-degree-of-freedom Hamiltonian systems given by time- 
independent Hamiltonians (one hrst integral reduces the dimension of the phase space to 3 so that the 
dimension of the Poincare hypersurface is equal to 2) and a wealth of information on the dynamical 
behaviour of the original system can be often obtained from this visualization. However, we consider 
Hamiltonian systems with N degrees of freedom with A > 3 in this paper, and, for this reason, there 
is no reasonable way how to visualize the Poincare hypersurfaces (their dimension is equal to 2N — 2), 
and hence it would require much more ehort to obtain any information on the behavior of the original 
system, which is beyond the scope of the present article. 

5. Concluding remarks 

To sum up, we have proved that, for an arbitrary even number k > 4, the system ([T]) with Hamiltonian 
@ does not admit an additional hrst integral meromorphic in a complex neighbourhood of the linear 
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manifold Qm = dmlN: Pm = o^mPN, m = 1,..., — 1, am = (—for m odd, am = 0 for m even, 

for all > 3 odd such that A^ = 1 mod 6 or A^ = 3 mod 6, and that the same system admits at 
most four additional first integrals meromorphic near the linear manifold mentioned above for A^ > 11 
odd such that N = 5 mod 6 (the case N = 5 with /c > 4 was settled by Maciejewski, Przybylska and 
Yoshida in m)- Moreover, computer-aided numerical experiments suggest that even in this case (i.e. 
N = 5 mod Q, N > 11), the system has no additional meromorphic first integral for “most” of the 
values of fc > 4 even. To give a more decisive answer, we have to solve two different problems: first, 
find all A: > 4 even, such that 5k'^ — 8k + 4 and 13A:^ — 16k + 4, are perfect squares, respectively, and 
second, even if for inhnitely many /c’s the corresponding values of the expressions — 8k + 4 and 
13A:^ — 16k -|- 4 are perfect squares, it could well happen that the original system does not admit an 
additional meromorphic first integral as well, the poly-Painleve test being simply indecisive here, so 
that still different tests are needed. Both of these last two problems deserve to be studied separately. 

Last but not least, a natural question arises: what happens if we consider different values of N 
and k than those considered above. It turns out that for k = 0,1,2 (with N arbitrary) the corre¬ 
sponding system is linear and therefore integrable in any reasonable sense. For A^ > 0 even (with an 
arbitrary fc > 3) we weren’t able to find any nonzero point a in a closed form for which the proce¬ 
dure would work. For N > 5 odd that is congruent to 2 mod 3 (with k > 3 odd) we found a point 
a = {a, —a, 0, a, —a, 0,..., a, —a) , where a = {k{l — However, in this case we weren’t 

able to say too much about the eigenvalues of the associated matrix Mp in general. We made some 
numerical experiments for N = 5 and k >3 odd and it seems that for k>21 all the four numbers Vm 
could be irrational so that the system ([2]) seems to possess no additional first integral meromorphic in a 
complex neighbourhood of the linear manifold Qm = Pm = ^Pn, m = 1,..., N — 1 in this case. 

In other cases of N and k remaining, we weren’t again able to find the point a. 

The results used in this paper rest essentially on the so called poly-Painleve method which was 
proposed at the beginning of 90’s by M. Kruskal. This method is based on asymptotic expansions 
of the unknown solutions, an idea which is close to the so called Painleve a-method. It seems that 
Kruskal himself in fact conjectured that if dense branching of a truncated asymptotic series to a certain 
order occurs, the actual solutions also have a dense branching and therefore there are no continuous 
first integrals of the corresponding equation. Although this conjecture is still not proved in its full 
generality, there are some partial results where this method is rigorously justified 0.0. a. 0. one of 
which was used in our proof. 
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